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O 

o 

< 

O 

cn 



Katsunori KuBO* and Yoshio Kuramoto 

Department of Physics, Tohoku University, Sendai 980-8578 
(Received April 21, 2003) 

Possible order parameters of the phase IV in CexLai-^^e are discussed with special attention 
to the lattice distortion recently observed. A Fsu-type octupole order with finite wave number 
is proposed as the origin of the distortion along the [111] direction. The Fg crystalline electric 
field (CEF) level splits into three levels by a mean field with the Fs^ symmetry. The ground 
and highest singlets have the same quadrupole moment, while the intermediate doublet has an 
opposite sign. It is shown that any coUinear order of Fsu-type octupole moment accompanies 
the Fsg-type ferro- quadrupole order, and the coupling of the quadrupole moment with the lattice 
induces the distortion. The cusp in the magnetization at the phase transition is reproduced, 
but the internal magnetic field due to the octupole moment is smaller than the observed one 
by an order of magnitude. 
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The multipolar orderings in cubic rare-earth hexa- 
borides have been studied extensively. Especially, CeBg 
is a typical material which shows multipolar orderings. 
In CeBg, there are three phases: the paramagnetic phase 
(called the phase I), the antiferro-quadrupole phase (II) 
and the antiferromagnetic phase (III). In Ce2,Lai_j;Bg, 
another phase, so called phase IV, was found at x ~ 
0.75."'^^ This phase has attracted much attention, but 
the order parameter has not yet been established. In the 
phase IV, a large softening of the elastic constant C44 
was observed. Furthermore the magnetic susceptibil- 
ity shows a cusp at the transition temperature Ti-iv from 
the phase I, and the magnetization is almost isotropic at 
ambient pressure. 

Recently, a small lattice distortion along the [111] di- 
rection was observed in the phase IV.^'^^ It is proba- 
ble that the softening of C44 is associated with this dis- 
tortion. It may be tempting to ascribe the distortion 
to the ferro Tsg-type quadrupole order. However, the 
quadrupole moment is not necessarily the primary order 
parameter. Large change of the internal field at Ti.rv, 
as probed by NMR^^ and fiSKj^ suggests strongly that 
the time reversal symmetry is broken in the phase IV, 
which is incompatible with a pure quadrupole order. Be- 
sides, neutron diffraction experiment found no magnetic 
reflection in the phase IV. Thus dipole moments are 
unlikely to be the primary order parameter either, al- 
though the time reversal symmetry is broken. Therefore, 
octupole moments, which break the time reversal sym- 
metry, become a candidate for the order parameter in 
phase IV. ^' ^^"^ Since the cubic symmetry is broken, and 
since the anisotropy develops in the magnetization under 
uniaxial pressure, the order parameter should have 
an anisotropic nature. Thus the Tsu-type, among all oc- 
tupole moments, is the most plausible candidate for the 
order parameter in the phase IV. 

Kusunose and Kuramoto^*^-* have already pointed out 
using the Ginzburg-Landau (GL) theory that the T^u 
octupole order with finite wave number should accom- 
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pany a ferro-quadrupole moment, and have suggested a 
possible lattice distortion. However, evaluation of the 
magnitude of the distortion is beyond their GL theory. 
In this paper, we explore in much greater detail the con- 
sequence of the Tsu octupole order by the mean field 
theory, and propose that the lattice distortion is due to 
the order of the Tsu octupole moment. 

The CEF ground state of Ce^+ (J = 5/2) in 
CcajLai-j^Bg is the Tg quartet. ^'^"^^^ The excited level 
Tj lies about 500K from the Fg level and is neglected. 
The Fg states are represented in terms of eigenstates of 
Jr as 



|+T) = ^5/61+5/2) 
|-T> = l + l/2), 



3/2), 



(1) 
(2) 



where + and — denote orbital indices, and their Kramers 
partners |+ i), ]— i) are obtained by reversing the signs 
of Jz in eqs.(l) and (2), respectively. Within the Fg quar- 
tet, the number of independent multipolar moments is 
15, and active octupole moments have either of F2u, F4„ 
or F5„-type symmetry. The F4„-type octupole mo- 
ments accompany dipole moments, and are unlikely 
to be the order parameter in the phase IV. The F2„-type 
octupole moment, as proposed by ref.9, does not accom- 
pany quadrupole moments. Therefore it seems difficult 
to explain the lattice distortion in the phase IV by the 
F2„ octupole orderings. 

Let us introduce pseudospins cr and t to describe the 
Fg quartet: 



(3) 



a^|r^t) = +|r^T), <j'\r' i) = -\r' i) , (4) 
and the transverse components which flip the pseu- 
dospins. The F5„ octupole, T^g and T^g quadrupole mo- 
ments are given with the notation (a,/3, 7) = {x,y,z), 
(y,z,x) or {z,x,y) by 



p - J2JJ/(2V3) = (C+a^Ca^rV^), 



(5) 
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Fig. 1. The level scheme in the Fs^ interaction (right) and in the 
Fsg interaction (left). 



(2J|- Jj- j2)/y3=(8/V3K, 
J,2-J,2 = (8/x/3K, 
0;3^ = = (2/\/3)r%«, 



(6) 
(7) 
(8) 



where bars on the products represent symmetrization, 
e.g., JxJl = (JxJy + JyJxJy + JyJx)/i, and we 
have introduced the notation = — (r^ ± •\/3t^)/2. 
We find that the 'easy axis' of T^u octupole moment 
is along the [111] direction and equivalent ones, i.e., 
(|(T5«)|, |(T5«)|, |(T5«)|) II (1, 1, 1) provided the intersite 
interaction is isotropic. Then we consider the following 
single-site Hamiltonian: 



n_J A (rpiiu I rp^u i rp^u\ 

'Isingle-site — -^5ul,-'x ' -'■ y ' -'■ z ) 



(9) 



where A^u denotes the octupolar mean field. The energy 
levels of this Hamiltonian are shown in the right part of 
Fig. 1. The Fg level splits into three levels. Not only the 
time reversal symmetry, but also the cubic symmetry 
are broken. The multipolar moments in the ground and 
highest states are (T^") = (T^") = (T^") ^ Ty^' 
{Oyz) = (Ozx) = {Oxy) = —2/3 and the others are zero. 
Thus, ferro, antiferro and other coUinearly ordered states 
with ((T5«), (T^''), (T^")) II (1, 1, 1) have a homogeneous 
r^g moment and the crystal distorts along [111]. 

Obviously a Fsg-type ferro-quadrupole interaction can 
also lead to {Oyz) = {Ozx) = {Oxy) 7^ 0. In this case 
the energy level splits into two levels, and the ground 
state is two-fold degenerate as shown in the left part of 
Fig. 1. Experimentally, the entropy changes much at 
the transition from the phase I to IV, and little from 
the phase IV to III in Ceo.75Lao.25B6.^' Thus it is 
likely that the degeneracy of the each /-electron state is 
already lifted in the phase IV. This situation also makes 
it hard for the T^g quadrupoles to be the order parameter 
in the phase IV. 

From the above discussion, we find that the order pa- 
rameter which breaks the time reversal symmetry and 
induces quadrupole moments, but which accompanies no 
dipole moment is only the octupole moment. How- 
ever, we cannot determine the periodicity of the F5„ or- 
dered state in the phase IV only from the above con- 
sideration, because any coUinearly ordered state with 
((r5-),(T5«),(r5")) II (1,1,1) has a uniform T^g mo- 
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Fig. 2. Temperature dependence of the antiferro-octupole moment 

(r5«> (= (tS") = (T|«)). 




Fig. 3. Temperature dependence of the ferro-quadrupole moment 

{Oy^) (= (O^x) = (Ox„». 



ment. To carry out the mean field theory explicitly, we 
consider a G-type antiferro-octupole order as the sim- 
plest example of F^u orders. The importance of 
nearest-neighbor interaction in causing the change from 
the phase III to III', even with weak magnetic field, was 
pointed out in ref.lO, and we consider only this nearest- 
neighbor interaction. We take the following model: 



n 



E 



E 



a i a j ' 



(10) 



where denotes a nearest-neighbor pair. We study 
this Hamiltonian by the mean field theory, and choose 
the value of so as to reproduce the transition tem- 
perature Tmv in Ceo.75Lao.25B6, i.e., T^u = 67^" = 1.7K. 

The temperature dependence of (T^") and {Oyz) 
are shown in Figs. 2 and 3. The solution ob- 
tained has the antiferro-octupole order with (T^") = 
{Ty^) = {T^"), accompanying ferro-quadrupole moment 
{Oyz) = {Ozx) = {Oxy)- We also find other equivalent 
solutions: ((T^"), (T|«), (T^"), (O^^), (0„^)) = 

{±B,±B,tB,+C.+C,-C) where quantities B and C 
depend on temperature. 

In Fig. 4, we show the magnetization in magnetic field 
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Fig. 4. Temperature dependence of the magnetization in magnetic 
field = 0.2T along various directions. 



H = 0.2T along three high-symmetry directions. The 
magnetization has a cusp at T5„, which is consistent 
with experimental observation.''^ We have also exam- 
ined the case where the pure quadrupole moment is 
the order parameter, and found that the magnetization 
changes little at the quadrupole transition temperature. 
The magnetization below Tsu is anisotropic. One should 
not, however, take this anisotropy seriously since other 
interactions such as quadrupole and dipole interactions 
also influence the anisotropy. In fact the easy axis in the 
mean field theory is different from that observed in the 
phase IV under uniaxial pressure and in the phase III. 

We now consider lattice distortion in the antiferro- 
octupole ordered state. The elastic energy associated 
with the moments is given per unit volume by 



E = 




(11) 

where eafi is the strain tensor, grr, is the magneto-elastic 
coupling constant, and C^^l^ is the (bare) elastic constant. 
The sum runs over Ce sites i in the unit volume, 
minimizing the elastic energy, we obtain 



By 



{Oa/3 i)- 



(12) 



with {Oy, ,) = (O^^ ,) _= {O^y ,) = -2/3. Note that 
along directions [111], [HI], and [111] the lattice should 
expand by 



M/l = 2\ey,\/'A = 3.1 X 10" 



(15) 



We use the following experimental values for Ceo.75Lao.25B6 
Iffrs I = 155K, Cf^ ~ 8.2 x lOiierg/cm^,^) and the lattice 
constant a = 4.13A. At absolute zero, the magnitude of 
quadrupole moments is given by \{Oyz i)| = \{Ozx i)\ = 
\{Oxy i)\ = 2/3, and we obtain 



4.6 X 10" 



(13) 



In order to account for the observed lattice contraction 
along [111], we consider the two possibilities: (i) gr^ < 0, 
and (ii) > 0. In the case (i), a slight stress accom- 
panying the measurement breaks the equivalence of four 
octupole domains, and choose a domain for which the 
contraction becomes the maximum along [111]. Namely 
we have 



Al/l = 2e 



-9.3 X 10" 



(14) 



On the other hand, in the case (ii), a positive stress along 
[111] may favor a domain which contracts along this di- 
rection. Then the contraction is given by eq.(15) with the 
minus sign. In the single domain with ~{Oyz) = (Ozx) = 
{Oxy) for example, an expansion along [111] should be 
present, and its magnitude is three times larger than the 
contraction along [111]. 

A magnetic field induces antiferromagnetic moments 
which tend to be perpendicular to H . Hence with H \\ 
[111] without external stress, a state where one of {Oai3)'s 
has a sign diff'erent from the others is stabilized. Exper- 
imentally the contraction is enhanced under H \\ [111], 
which alone favors the case (ii) in our model. In the 
actual system, diagonal components becomes posi- 
tive in the phase IV. This comes from the volume strain 
which is not included in our model. According to ex- 
periment,^^ the shear strain with the assumption of the 
trigonal symmetry aroimd [111], i.e., the case (i), is de- 
rived as ea/3 = — 4x 10~^ without magnetic field at 1.3K. 
Assuming that the phase IV is stable down to zero tem- 
perature, extrapolate to (—6 ^ —10) x 10~^ at abso- 
lute zero. The absolute value is by an order of magnitude 
smaller than our estimate in eq.(13). Then the observed 
reduction of |(0a/3)| should come from quantum fluctua- 
tions neglected here. Another possibility is that the case 
(ii) is realized. In this case three times larger |e„/3|, i.e., 
(2 ^ 3) X 10~^, is obtained from the same experimental 
result.^) Our model in the case (ii) predicts a larger lat- 
tice contraction along the [110] direction with H || [110] 
than that along the [111] direction with H \\ [111]. 

For comparison, we estimate in the same way the mag- 
nitude of the lattice distortion of NdBe in the O2 ferro- 
quadrupole ordered state by using the experimental val- 
ues: Iffral = 220K, (Cll-Cl2)/2~20xl0"erg/cm^l8) 
a = 4. 12 A, and the Lea-Leask-Wolf parameter sjllw = 
-0.82.19) We assume that (cj"^ - cf2')/2 is almost 
the same as the observed (Cn — Ci2)/2. In a magnetic 
field H II [100], a state with the antiferromagnetic mo- 
ment perpendicular to [100] is stabilized. For a domain 
where the magnetic moment is along [001], we obtain 
|AZ|/Z = 2.8 x 10~^ along H. This value compares fa- 
vorably with the experimental one Al/l ~ —2 x 10~^ at 
2K and H = 2.IT.20) The distortion in NdBg is by an 
order of magnitude larger than that in eqs.(14) and (15). 
The reason is that the quadrupole moment 1(0") | = 4.5 
in NdBg is much larger than the corresponding value 
|(Oj/z)| = 2/3 in the octupolar state of CcxLai-xBg, 
while the magneto-elastic coupling constants are of the 
same order. 

We now discuss the internal magnetic field associated 

with the octupole order. The /xSR time spectra in the 
phase IV consist of a Gaussian component and an ex- 
ponential component.^) The observation of a Gaussian 
relaxation indicates that internal fields are randomly dis- 
tributed, and/or fluctuating. The internal field deduced 
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from the Gaussian relaxation is the order of O.IT. We 
discuss whether the octupole moment can be the ori- 
gin of the Gaussian relaxation. In /xSR measurement, 
/x"*" locates at (a/2,0,0) and equivalent sites^^"^"^) with 
a Ce ion chosen as the origin. As a reference the in- 
ternal field from a Bohr magneton fis is estimated to 
be i/dipoic = HB/{a/'^f ~ O.IT, with a/2 ~ 2A. The 
internal field from an octupole with the size r is esti- 
mated to be i?octupoie = Mb^^/('2/2)^ ~ O.OIT, with 
r ^ Ob = 0.53A. For a more accurate estimate, wc con- 
sider the multipolc expansion of the vector potential from 
local electrons as given by^^' 

^w = ET^"^'^'^O^™H^,0))Mr, (16) 

where I is the orbital angular momentum operator, 
Cm {0, (j)) is ^ All/ {2k + 1) times the spherical harmon- 
ics Yfem(^^, 0), and M™ is the magnetic multipole mo- 
ment. The multipole moment is determined by the wave 
function il)i{r) of the i-th electron, the orbital and spin 
angular momentum operators Z, and s,. Namely we have 

' (17) 

Eq.(17) is evaluated with use of the operator equivalents 
method. For our purpose it is sufficient to consider 
only octupole moments. Then we obtain through 
calculation of the reduced matrix element of the third 
rank tensor. The result for one electron states with J = 
5/2, L = 3, 5 = 1/2 is given by 

M^ = ~^,n{r'){J(^^), (18) 

(3) 

where the third-rank tensor operators Jm are defined in 
ref. 16. 

Freeman and Dcsclaux obtained the estimate (r^) = 
1.298 in atomic unit by a relativistic Dirac-Fock calcula- 
tion.^^) By using this value, we obtain about 40G as the 
magnitude of internal field at (a/2,0,0) from a Cc ion. 
This value is by an order of magnitude smaller than that 
derived by the /itSR measurement. Thus the static oc- 
tupolar moment alone cannot account for the Gaussian 
relaxation of the /uSR spectra. In a future work, we plan 
to study in more detail fluctuations in the octupole or- 
dered state. 

We mention that Paixao et al. have proposed for 
the ordered state of Np02 that a triple-g Fs^-type oc- 
tupole ordering is realized and a triple-q Fs^ quadrupole 
moment is induced. We estimate the internal mag- 
netic field in Np02 to be about 80G at 8K from the ob- 
served muon spin precessing frequency of TMHz.^^^ The 
internal field is of the same order as our estimate for 
Ce^Lai-xBe. 

In the antiferro-octupolar phase, antifcrromagnctism 
should be induced by uniaxial stress.^' In the following 
we estimate the magnitude of the induced moment. The 
uniaxial pressure p along the [001] direction accompanies 



the F3 strain 

(Cii-Ci2)(2e,^-e^^-ej,j,) = 2p. (19) 

We have solved the mean field equation with the finite 
F3 strain and the corresponding quadrupolc-strain inter- 
action g^^ . The antiferromagnetic moment is induced in 
the xj/-plane since (T^") 7^ together with ( - ) ^ 
gives {Jx) 7^ 0. The direction of the antiferromagnetic 
moment is along the [110] or [110]. With experimen- 
tal values: (Cn - Ci2)/2 ~ 20.4 x 10i^(erg/cm3) and 
\gT:,\ = 120K,2) the mean-field solution for (Tj^") be- 
comes zero for p > 0.9GPa. This implies the octupole 
and dipole moments both lying in the a;y-plane. With 
p — IGPa, the magnetic moment is estimated to be 
0.88/xb (0.82/xb), if gv^ is positive (negative). This value 
should actually be reduced by quantum fluctuations. 
However, it is likely that the magnitude remains within 
experimental access. 

To summarize, we have proposed that the F5„ octupole 
moment is a plausible candidate for the order parame- 
ter of the phase IV. This ordered state is consistent with 
the lattice distortion along the [111] direction,*'^) the 
broken time reversal symmetry,^' no dipole moment 
found, ^) and the cusp in the magnetization at Tj.iy."^) 
However, the internal field estimated by our mean field 
theory is much smaller than that suggested by the /iSR 
experiment.^-' We recall that the NMR spectra become 
broad in the phase IV, and the spectra cannot be ex- 
plained by the static and staggered octupole moments 
either. Therefore clarifying dynamical aspects and iden- 
tifying the periodicity in the phase IV are challenging 
open problems to be addressed in the near future. 
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